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ABSTRACT. The optimal dividend problem by De Finetti (1957) has been recently generalized to the spec- 
trally negative Levy model where the implementation of optimal strategies draws upon the computation of 
scale functions and their derivatives. This paper proposes a phase-type fitting approximation of the optimal 
strategy. We consider spectrally negative Levy processes with phase-type jumps as well as meromorphic 
Levy processes (Kuznetsov et al., 2010a), and use their scale functions to approximate the scale function 
for a general spectrally negative Levy process. We obtain analytically the convergence results and illustrate 
numerically the effectiveness of the approximation methods using examples with the spectrally negative 
Levy process with i.i.d. Weibull-distributed jumps, the /3-family and CGMY process. 
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1. Introduction 

In the optimal dividend problem, an insurance company wants to maximize the cumulative amount of dividends 
paid out to its beneficiaries until the time of ruin. The original problem by De Finetti (1957) considered a discrete- 
time case by modeling the surplus by a random walk. It has been later extended to the continuous-time diffusion 
model (Jeanblanc and Shiryaev, 1995; Asmussen and Taksar, 1997; Gerber and Shiu, 2004) and to the Cramer- 
Lundberg model (Gerber, 1969; Azcue and Muler, 2005). It was recently generalized to the spectrally negative 
Levy model by Avram et al. (2007) and Kyprianou and Palmowski (2007) where the surplus is a general Levy 
process with only negative jumps. 
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The implementation of the optimal strategies in the spectrally negative Levy model draws upon the computation 
of the so-called scale function. Avram et al. (2007) first obtained the expected value under a barrier strategy in 
terms of the scale function. Loeffen (2008) then showed that a barrier strategy is indeed optimal under a suitable 
condition. The scale function also plays great roles in its extensions with transaction costs and additional terminal 
values at ruin. See Loeffen (2009a,b) and Loeffen and Renaud (2010). 

Despite these advances, a major obstacle still remains in putting in practice the above-mentioned results because 
scale functions are in general known only up to their Laplace transforms. The implementation is even harder for 
the optimal dividend problem because it requires the first derivative of the scale function. One can in principle 
approximate the scale function and its derivative via the numerical Laplace inversion as discussed in Surya (2008). 
However, the method is only heuristic and one cannot determine the accuracy of the approximation. Moreover, the 
derivative of the scale function tends to explode in the neighborhood of zero, and therefore the error is expected 
to be large near zero. For these reasons, there is a clear need of new approximation methods that work for any 
spectrally negative Levy process. 

As a new tool to overcome these problems, this paper proposes a phase-type fitting approach to approximate the 
scale functions as well as the solutions to the optimal dividend problem and its extensions. We obtain the scale 
functions of spectrally negative Levy processes with phase-type jumps and those in the M-class (meromorphic 
Levy processes), and apply these for a general spectrally negative Levy model. 

We first consider the class of spectrally negative Levy processes with phase-type jumps. Consider a continuous- 
time Markov chain with some initial distribution and state space consisting of a single absorbing state and a finite 
number of transient states. The phase-type distribution is the distribution of the time to absorption. It is known 
that the class of phase-type distributions is dense in the class of all positive-valued distributions. We obtain the 
scale functions for these processes and show that they can approximate the scale function of a general spectrally 
negative Levy process arbitrarily closely. 

The phase-type fitting approach has mainly three advantages. First, thanks to the smoothness and monotonicity 
properties of scale functions as proved by, for example, Chan et al. (2009) and Loeffen (2008), the approximation 
can be applied also to its derivative. Second, the Laplace transform of the phase-type distribution has an explicit 
expression and hence can avoid the error caused while approximating the Laplace transform for a general jump 
distribution. This type of errors tends to occur in other approximation methods such as Surya (2008). Third, the 
phase-type fitting approach enjoys a variety of fitting algorithms. See, for example, Asmussen (1996), Bladt et al. 
(2003) and Feldmann and Whitt (1998). The fitting can be applied also to empirical data and this is another major 
advantage. In the first half of our numerical results, we use the results by Feldmann and Whitt (1998) and consider 
the case with a Brownian motion plus a compound Poisson with Weibull-distributed jumps. 
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We also consider meromorphic Levy processes (Kuznetsov et al., 2010a), which generalizes a number of 
recently-discovered Levy processes such as Lamperti-stable processes (Caballero and Chaumont, 2006; Chau- 
mont et al., 2009), hypergeometric processes (Kyprianou, 2010; Kuznetsov et al., 2010b) and processes in the (3- 
and ©-families (Kuznetsov, 2009a,b). The spectrally negative versions of these processes commonly have Levy 
measures in the form 



which can be seen as an extension of the hyperexponential density and as a "discrete version" of the completely 
monotone density. 

The reason we consider the Levy measure above is that it can approximate efficiently any Levy measure with 
a completely monotone density. Furthermore, if the Levy measure is completely monotone, the barrier-strategy is 
guaranteed to be optimal for the dividend problem (Loeffen, 2008). The class of Levy measures with completely 
monotone densities is rich. It enables us to model compound-Poisson-type jumps with long-tailed distributions 
such as the Pareto, Weibull and gamma distributions; see Feldmann and Whitt (1998). It further allows us to 
construct many of recently-introduced Levy processes such as variance gamma processes (Madan and Milne, 
1991; Madan et al., 1998), CGMY processes (Carr et al., 2002), generalized hyperbolic processes (Eberlein et al., 
1998) and normal inverse Gaussian processes (Barndorff-Nielsen, 1998). These processes can be approximated 
efficiently by meromorphic Levy processes. 

Given a Levy measure in the form (1.1), the corresponding scale function and its derivative can be expressed 
explicitly as infinite sums of exponential functions. We first show that these can be approximated by finite sums 
with some analytical error bounds, concluding that the error bounds for the solutions to the optimal dividend 
problem can also be obtained. We then show numerically the effectiveness of the approximation procedure using 
examples with /3-processes (Kuznetsov, 2009a) and CGMY processes. We obtain bounds on scale functions and 
solutions for the former and use them to approximate for the latter. 

It should be emphasized here that the approximation procedure discussed in this paper can be applied outside 
the class of optimal dividend problems. The spectrally negative Levy model has been recently introduced widely 
and the scale function plays a great role in characterizing the solutions. We refer the reader to Avram et al. (2004) 
and Alili and Kyprianou (2005) for derivative pricing, Kyprianou and Surya (2007) for optimal capital structure, 
Baurdoux and Kyprianou (2008, 2009) for stochastic games. For a comprehensive account, see Kyprianou (2006). 

The rest of the paper is organized as follows. Section 2 summarizes the results on the classical dividend problem 
and its extensions. In Section 3, we obtain the scale functions for spectrally negative Levy processes with phase- 
type jumps and show that it can approximate the scale function of a general spectrally negative Levy process 



(1.1) 
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arbitrarily closely, 
and lower bounds, 
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We then obtain in Section 4 the scale functions for those in the M-class including its upper 
We conclude this paper by giving numerical results in Section 5. All proofs are given in the 



2. Optimal Solutions to Dividend Problem via Scale Functions 

This section reviews the classical dividend problem and its extensions focusing on the spectrally negative Levy 
model. Let (ft, J 7 , P) be a probability space hosting a spectrally negative Levy process X = {X t ; t > 0} with its 
Laplace exponent 

1 f° 

(2.1) :=logE[e sXl ] = /xs + -aV + / (e sz - I - szl {z> ^ 1} )u(dz), sGC 

where v is a Levy measure with the support (— oo, 0) that satisfies the integrability condition J,^ ^ (1 Az 2 )v(dz) < 
oo. Moreover, let F x be the conditional probability under which Xq = x (also let P = P°), and F := {Ft : t > 0} 
be the filtration generated by X. The process X is called the Levy insurance risk process (or the risk process in 
short), and models the surplus of an insurance company before dividends are deducted. 

The classical dividend problem is a control problem where the cumulative amount of dividends prior to ruin is 
maximized. A (dividend) strategy ir := {LJ; t > 0} is given by a nondecreasing, left-continuous and Y-adapted 
process starting at zero. Corresponding to every strategy it, the remaining amount of surplus after dividends are 
deducted is given by U n = {Uf : t > 0} where 

U?:=X t -I%, t>0, 

and its ruin time is the first time it goes below zero: 

a" = inf {t > : U? < 0} . 

A lump-sum payment must be smaller than the available fund and hence it is required that 

(2.2) L? + - < U?, t < a w a.s. 

Let n be the set of all admissible strategies satisfying (2.2). The problem concerns the expected sum of total 
discounted dividends until ruin 



vJx) := E* 



e- qt &Ll 



o 



i s n, 



and wants to obtain an admissible strategy that maximizes it. Hence the classical dividend problem is written as 

(2.3) v{x) := sup 1^(0:). 

vren 
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2.1. Scale functions. The solutions to the dividend problems can be written in terms of the scale function. Here, 
we describe the scale function and summarize its properties that will be used in this paper. 

Fix q > and any spectrally negative Levy process with its Laplace exponent tp as defined in (2.1). The scale 
function W {q] : is a function whose Laplace transform is given by 

1 



(2.4) 
where 
(2.5) 



oo 

e- sx W iq \x)dx = — — 



S > Cq 



C q := sup{s > : ip(s) = q}, q> 0. 



We assume W<d(x) = Oon (-oo,0). 

Let us define the first down- and up-crossing times, respectively, by 



(2.6) 



r a : = inf {t > : X t < a} and r+ := inf {t > : X t > b} , 



for every < a < x < b. Then we have 

W^(x) 



(2.7) E 
where 
(2.8) 



e- qT b 1 



W<ro} 



and E x 



{r+>r } 



zte\x)-z®(b) 



W^\x) 
W(i)(b) 



Z^ q \x) := 1+q [ W iq) {y)dy, x£K. 
Jo 



Here, we disregard the case when X is a negative subordinator (or decreasing a.s.). 
We also consider a version of the scale function Wq that satisfies 



(2.9) 

with its Laplace transform 
(2.10) 



W (q) {x) = e c " x W Cq (x), i£ 



OO 1 



s > 0. 



Suppose P c , for any given c > 0, is the probability measure defined by the Esscher transform 



dP 



see page 78 of Kyprianou (2006). Then under is analogous to under P. Furthermore, it is known 
that W^ q is monotonically increasing and 



(2.11) 



W Cq {x) SWiQT 1 asx 



oo, 
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which also implies that the scale function increases exponentially in x; 

C x 

(2.12) W iq) (x) ~ 4ttt\ asx^oo. 

Due to the fact that Wq does not explode for large x as opposed to W^ g \ it is often convenient to deal with Wq 
and convert it to W^ 1 ' using (2.9), especially when numerical computations are involved; see Surya (2008). 
Recall that a spectrally negative Levy process has paths of bounded variation if and only if 



<7 = and / (1 A |x|)^(dx) < oo; 

J(-oo.O) 



'(-oo,0) 

see, for example, Kyprianou (2006), Lemma 2.12. In this case, we can rewrite the Laplace exponent (2.1) by 



ip{s) = lis + \a 2 s 2 + f (e sz - l)v(dz), s G C, 



with 



fO 

p. := fi — I zu(dz). 



-l 

Regarding the smoothness of the scale function, it has been shown by Chan et al. (2009) that if a Levy process 
has a Gaussian component (a > 0), we have W^ q ' G C 2 (0, oo). When it does not have a Gaussian component 
and if its jump distribution has no atoms, we have G C 1 (0, oo). In particular, a stronger result holds for 

the completely monotone jump case. Recall that a density function / is called completely monotone if all the 
derivatives exist and, for every n > 1, 

(-l) n / (n) (*)>0, x>0, 

where denotes the n th derivative of /. 

Lemma 2.1 (Loeffen (2008)). If the Levy measure has a completely monotone density, is again completely 
monotone. 

Finally, the behavior in the neighborhood of zero is given as follows. See Lemmas 4.3 and 4.4 of Kyprianou 
and Surya (2007). 

Lemma 2.2. For every q > 0, we have 

'hi ^>o 



i 0, unbounded variation ■ 
W^iO) = { } and W {q) (0+) = < 

i , bounded variation 



oo, (7 = and u(—oo, 0) = oo 

q+v ^~£°'^ , compound Poisson 
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2.2. Solutions in terms of scale functions. In a spectrally negative Levy model, the expected value under the 
barrier strategy can be expressed in terms of the scale function as shown by Avram et al. (2007). A barrier strategy 
at level a > is denoted by ir a := {Lf ; t < a a } where 

L a t := sup(X s - a) V 0, t > 0. 

s<t 

We further let a a '■= inf {t > : Vj a < 0} denote the corresponding ruin time. 

Theorem 2.1 (Avram et al. (2007), (5.1)). For every a > 0, we have 

u a (x), x < a, 



(2.13) v Va (x) = 
where 

(2.14) u a {x) : 



x — a + u a {a), x > a, 
W (q \x) 



Furthermore, the barrier strategy attains optimality under a suitable condition. 

Theorem 2.2 (Loeffen (2008)). Suppose a* satisfies 

W^ q) \a) < W (q) \b), a* <a<b. 
Then the barrier strategy ir a * is an optimal strategy. 

In view of the above, an optimal barrier strategy exists, for example, when W' is convex, which holds whenever 
Levy measure has a completely monotone density as in Lemma 2.1. 

Extension with bail-out. A variant called the bail-out problem is discussed in Avram et al. (2007). Here the 
beneficiary of the dividends must inject capital to keep the risk process from going below zero. A strategy is now 
a pair 7f = {L w , i? 71 "} where L w is the cumulative amount of dividends as in the classical model and R w is a 
right-continuous process representing the cumulative amount of injected capital satisfying 



(2.15) / e~ qt dB? < oo 



a.s. 



o 



Assume that ip > 1 is the cost per unit injected capital, the problem is to maximize 

v w (x) = E 2 



e~ qt dLJ - ^ / e- qt dK" t 
ii Jo 



among all strategies that satisfy (2.2) and (2.15). In this model, the optimal strategy reduces to the double-barrier 
strategy that regulates the risk process Vt := X t — L\ + inside the interval [0, d*] where 

(2.16) d* :=inf {a > : G{a) := (<pZ {q \a) - 1)W W (a) - pqW^ia) 2 < o} . 
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The value function becomes Vd* (x) with 

x - d + v d (d), 

Extension with terminal values at ruin. Loeffen and Renaud (2010) recently considered the case with addi- 
tional terminal costs at ruin; the objective function is 



(2.17) 



Vd{x) 



qWW(d) 



, 0<X<d, 

x > d. 



maxE x 
vren 



{a *<oo) 



for some affine function P(y) := S + Ky. Its special case with constant terminal value was studied by Loeffen 
(2009a). With the assumption that the tail of the Levy measure is log-convex, the optimal solution is either the 
barrier strategy or the take-the-money-and-run strategy where the latter immediately pays out all the dividends and 
forces the ruin to occur immediately. For the former case, the optimal barrier level is given by 



b* := sup {b > : F(b) > F(x) for all x > 0} 



where for every x > 

(2.18) 

and 



F{x) :-- 



1 - A(x) 
W(«Y(x) 



A(x) := K (z^ix) - i)'{0+)W {q) {x)} + SqW iq \x). 
The value function is given by u&* (x) with 

S + f x A(y)dy+^f ) WM(x), < x < b, 



(2.19) 



v b (x) 



x — b + Vb(b), x > b. 

Extension with transaction costs. An extension allowing transaction costs is discussed in Loeffen (2009b) 
where the objective is to maximize 



Jo \ 0<s<t 



>0} 



Here 6 > is the unit transaction cost and the strategy is assumed to be defined by a pure jump process in the form 



0<s<t 



For this impulse control problem, the role of the barrier strategy in the classical model is now replaced by the so- 
called (ci, C2)-policy which is commonly known in inventory control. The (ci, C2)-policy brings the risk process 
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down to the level c\ whenever the risk process goes above the level c 2 - Loeffen (2009b) showed that (cf , c^-policy 
is optimal if it satisfies c*2 — c\ — 5 > 0, minimizes the function 

, , W^(c 2 )-WW(ci) 

g{ci,c 2 ) = - , 

C2 — Cl — 

and satisfies W^' (a) < W®' '(b) for every c!^ < a < b. The value function has the form (2.13) where the optimal 
threshold level is replaced with c|. 

3. Scale functions for spectrally negative Levy processes with phase-type jumps 

As we have seen in the last section, the implementation of the optimal strategies in the optimal dividend problem 
draws upon the computation of scale functions. This section obtains the scale function of the spectrally negative 
Levy process with phase-type jumps and shows that it can be used to approximate the scale function of any 
spectrally negative Levy process. 

3.1. Spectrally negative Levy processes with phase-type jumps. Consider a continuous-time Markov chain 
Y = {Y t ;t > 0} with finite state space {1, . . . , m} U {A} where 1, . . . , m are transient and A is absorbing. Its 
initial distribution is given by a simplex a = [a±, . . . , a m ] such that a- L = P {Yq = i} for every i = 1, . . . , m. The 
intensity matrix Q is partitioned into the m transient states and the absorbing state A, and is given by 

T t 


Here T is an m x m-matrix called the phase-type generator, and t = —Tl where 1 = [1, . . . , 1]'. A distribution is 
called phase-type with representation (m,ac,T) if it is the distribution of the absorption time to A in the Markov 
chain described above. It is known that T is non-singular and thus invertible; see Asmussen (1996). Its distribution 
and density functions are given, respectively, by 

F(z) = 1 - ae Tz l and f(z) = ae Tz t, z > 0. 

Let X = {X t ;t > 0} be a spectrally negative Levy process of the form 

N t 

(3.1) X t - X = nt + aB t -J2 Z n, 0<t<oo, 

n=l 

for some \i € K and a > 0. Here B = {B t ; t > 0} is a standard Brownian motion, N = {N t ; t > 0} is a Poisson 
process with arrival rate A, and Z = {Z n ;n = 1,2,...} is an i.i.d. sequence of phase-type distributed random 
variables with representation (m,a,T). These processes are assumed independent. Its Laplace exponent is then 

(3.2) ij)(s) =/is+ \ 2 s 2 + A (a(sl - TyH - l) , 



Q 



10 M. EGAMI AND K. YAMAZAKI 

which is analytic for every s G C except for the eigenvalues of T. 

Disregarding the case when A" is a negative subordinator, we consider the following two cases: 

Case 1: when a > (i.e. X has unbounded variation), 

Case 2: when a = and \i > (i.e. X is a compound Poisson process). 

Notice, in Case 2, that we can write Xt = Ut — J2n=i %n where Ut = x + pi is a (positive) subordinator. This 
implies that down-crossing of a threshold can occur only by jumps; see, for example, Chapter III of Bertoin (1996). 
On the other hand, in Case 1, down-crossing can occur also by creeping downward (by the diffusion components). 
Due to this difference, the form of the scale function differs as we shall see. 
Fix q > 0. Consider the Cramer-Lundberg equation 

(3.3) il>(s) = q, 

and define the set of (the absolute values of) negative roots and the set of poles: 

X q := {i : ip{-Ci, q ) = q and TZ{^ q ) > 0} , 

J q ■= \j ■ Tl T = and > o) • 

The elements in Z q and J q may not be distinct, and, in this case, we take each as many times as its multiplicity. By 
Lemma 1 of Asmussen et al. (2004), we have 

1.7,1 + 1, for Case 1, 
\Jq\-, for Case 2. 

In particular, if the representation is minimal (see Asmussen et al. (2004)), we have \ J q \ = m. 

Let K q be an independent exponential random variable with parameter q and denote the running maximum and 
minimum, respectively, by 

Xt = sup X s and X_ t = inf X 8 , t > 0. 

0<s<t 0<s<t 

The Wiener-Hopf factorization states that q/(q — ip(s)) = tp+(s)ip~ (s) for every seC such that 1Z(s) = 0, with 
the Wiener-Hopf factors 

(3.4) <p-(s):=E[eMsX Kg )} and E [exp(«X^)] 

that are analytic for s with 1Z(s) > and TZ(s) < 0, respectively. By Lemma 1 of Asmussen et al. (2004), we 
have, for every s such that TZ(s) > 0, 



(3-5) <p-{s) 
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from which we can obtain the distribution of X_ Kq by the Laplace inverse via partial fraction expansion. 

As in Remark 4 of Asmussen et al. (2004), let n denote the number of different roots in Z q and mj denote the 
multiplicity of a root for i = 1, . . . , n. Then we have 

where 



(3-6) 



n mi , £ 1 



A 



(mi — k)l ds mi ~ k 



Notice that this can be simplified significantly when all the roots in X q are distinct. 

3.2. Scale functions for spectrally negative Levy processes with phase-type jumps. Here we obtain the scale 
function. We focus on the case q > because the scale function when q = (and X drifts to infinity) can be 
derived by using W^°\x) = F x {X^ > 0} /ip'{0) and the ruin probability (19) of Asmussen et al. (2004) by 
taking q — > 0. Kyprianou and Palmowski (2007) briefly stated the scale function when q = and all the roots in 
X q are distinct. 

Before obtaining the scale function, we shall first represent the positive root ( q (2.5) in terms of the negative 
roots {£i j(3 ; i G I q }. Let us define 

n 



i=i 



(3.7) 

and by Lemma 2.2 

(3.8) 6 := -( q W (q \0) + W {qY (0+) 

Lemma 3.1. For every q > 0, we have 

k 



'2 , for Case 1 

■$ L + q ¥, for Case 2 

11 1 it* ' 



Qq 



We now obtain the version of the scale function (•). In the lemma below, W^ q (0) = W^ q \0) is either or 
- depending on if it is Case 1 or Case 2; see Lemma 2.2. 



Lemma 3.2. For every q > 0, we have 

> n vrii ft \ fc 



l,q 



Cq + 



fc-l 



1 _ e -(c«+£i,«> ^2 

j=0 



((Cq + ^iM j 



* i=l k=l 

Lemma 3.2 together with (2.12) and Lemmas 2.2 and 3.1 shows the following. 



J! 



x > 0. 
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Proposition 3.1. For every q > and x > 0, we have the following. 
(1) For Case 1, we have 



(k) ( ihq 



k-1 



Cq x g £,i,q x ^ ^ 



((c,+&, g )*) j ' 



(2) For Case 2, we have 



W^\x) 



1 A Cg . 9 + A\ ^ ,(fc) / 



+ 



^ V ^ ^ / ,=i fc =i 



V V ( SM 

Z^ i,9 lc, + 6,« 



i=o 



((C, + &, 9 )*) 3 



+ -e c?x . 



Recall that the solution to the dividend problem requires the derivative. The scale functions obtained above are 
infinitely differentiable. In particular, the first derivative becomes 

fc-2 



t sr sr A( k ) ( 

a 2 n Z^Z_^ i,q \ (- _l i. 



for Case 1 and 



1//- i \ \ n m i / £ \ k 

■ / Cg | g + a ^ ^2 sjr A {k) ' ^ 



Qq \ H H 



1=1 fc=l 



+ (A ,_ 1)! ^ ^ ^ 



A* 



for Case 2. 

When all the roots in X q are distinct, the scale functions above can be simplified and have nice properties as 
discussed in the following corollary. 



Corollary 3.1. If all the roots in X q are distinct, we have the fallowings. 
(1) The scale function can be simplified to 



W (q \x) 
W( q \x) 



£4 



mf e 



i,q 



1 ( Cq , 9 + A 



Qq V A* A* 2 / ^ 



C fM ) 




\C + (i,,J 





A* 



/or Case 1 and Case 2, respectively. 

(2) W (q) ' is convex. 

(3) is completely monotone. 
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This guarantees the optimality of the barrier-strategy in view of Theorem 2.2; the optimal strategy can be 
obtained by finding a unique a* such that W^ q > (a*) = 0. 

Example 3.1 (Hyperexponential Case). As an important example where all the roots in X q are distinct, we consider 
the case where Z has a hyperexponential distribution with density function 

m 

f(z) = Y,^Vie-^ z , z>0, 
i=i 

for some < 771 < • • • < r] m < 00. Its Laplace exponent (2.1) is then 

- m 

(3.9) =/is + -aV-AV aj — . 

2 ~[ Vi + s 

Notice in this case that — r/i, . . . , —r] m are the poles of the Laplace exponent. Furthermore, all the roots in X q are 
distinct and satisfy the following interlacing condition for every q > 0: 

(1) when a > 0, there are m + 1 roots — • • • , — £m+l,g such that 

< 6,9 <m < 6,9 <"■ <Vm< €m+l,q < 

(2) when a = and fi > 0, there are m roots —£i q, • • ■ , —£m,q such that 

< £l,q < m < < • • • < S,m,q < 7J m < OO. 

The class of hyperexponential distributions is important as it is dense in the class of all positive-valued distributions 
with completely monotone densities. 

3.3. Approximation of the scale function of a general spectrally negative Levy process. The scale function 
obtained in Proposition 3.1 can be used to approximate the scale function of a general spectrally negative Levy 
process. By Proposition 1 of Asmussen et al. (2004), there exists, for any spectrally negative Levy process X, a 
sequence of spectrally negative Levy processes with phase-type jumps converging to X in D[0, 00). This 

(n) 

is equivalent to saying that X\ —¥ X\ in distribution by Jacod and Shiryaev (2003), Corollary VII 3.6; see also 
Pistorius (2006). Suppose ip n (ip), C g , n (C?) an d /W^ q! n (W^ /W^ q ) are the Laplace exponent, the positive 
root (2.5) and the scale function of X^ (X), respectively. Because these processes are spectrally negative and 
ip is continuous, we have, by the continuity theorem, tp n ((3 + Q,n) — > + C?) f° r every ft > 0. Now in 
view of (2.10), the convergence of the scale function holds by the continuity theorem; see Feller (1971), Theorem 
2a, XIII. 1. More precisely, we have Jj W^ qin (y)dy ->■ fj W^ q (y)dy and Jj wjf\y)dy — > J 7 W^ q \y)dy for any 
interval /. 
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The smoothness and monotonicity properties of the scale function can be additionally used to obtain stronger 
results. The scale functions in Proposition 3. 1 are in C°° (0, oo). In addition, when all the roots of X q are different, 
its first derivative is completely monotone as discussed in Corollary 3.1. 

If the target scale function is in C 1 (0, oo) (which holds whenever the jump distribution has no atoms), noting 
that Wq{x) < (^'(C?)) -1 f° r every x and hence e~^ x W^ q (x) vanishes in the limit for any j3 > 0, we have by 
(2.10) 



oo 

e 



Because W£ (x) is nonnegative and F(x) := Jq (y)dy / ((tp 1 (Cj)) 1 — Wc 9 (0)) * s a probability distribution, 

W Cq , n (x) ^> W Cq (x) and W${x) ^> W^(x), x>0. 

Furthermore, suppose that it is in C 2 (0, oo) (which holds, for example, when a > by Chan et al. (2009)), 
(0+) < oo (i.e., a > or II(— oo,0) < oo) and W" (x) < for every x > 0, because W^(x) ^2$. 0, we 
have F(x) := (Wi (0+)) _1 f£ \ W'J (y)\dy is a probability distribution and 



f 

Jo 



e 



-^F(dx) = (Wi(0+))- 1 



+ pw <q (o) + wi (0+) 



p > 0. 



+ d) 

Therefore, noting that ^'(i) = ( q W^(x) + e^^Ji) and assuming that the convergent sequence , n (x) 
has the same property, we can obtain by the continuity theorem 

Wl q>n {x) ^W' Cq {x) and W$' (x) ^ W iqY (x), x > 0. 

The negativity of W^' holds, for example, for the completely monotone jump case because is completely 
monotone by Loeffen (2008). We can also choose the sequence n completely monotone in view of Corollary 
3.1 because approximation can be done via hyperexponential distributions. In fact, it also means that is 
C°°(0, oo) and the convergence of higher derivatives can be pursued. Even for a general jump distribution, the 
negativity of is a reasonable assumption in view of the numerical plots given by Surya (2008). 

This phase-type fitting approach complements the approach by Surya (2008) where scale functions are approx- 
imated by numerical Laplace inversion. A major disadvantage of using this inversion method is the fact that it 
requires the exact value of the right-hand side of (2.4). However, the Laplace transform of a jump distribution does 
not in general have an explicit closed-form expression. Surya (2008)'s approach, therefore, contains two types of 
errors: 1) the approximation error caused while computing ip and 2) the error caused while inverting the Laplace 
transform. On the other hand, the phase-type fitting approach only contains the phase-type fitting error thanks to 
the closed-form Laplace transform of the phase-type distribution. 
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The phase-type fitting approach enjoys a variety of fitting algorithms typically developed in queueing analysis. 
Well-known examples are the moment-matching approach (e.g. MEFIT and MEDA) and the maximum-likelihood 
approach (e.g. MLAPH and EMPHT), and a thorough study of pros and cons of each fitting techniques has been 
conducted in, for example, Horvath and Telek (2000) and Lang and Arthur (1996). The fitting can be applied also 
to empirical data and this is another major advantage over the Laplace inversion approach. 

4. Scale functions for Meromorphic Levy processes 

In this section, we consider another class of spectrally negative Levy processes called meromorphic Levy pro- 
cesses. We obtain their scale functions and use these as approximation tools for a general spectrally negative Levy 
process with a completely monotone Levy measure. Similarly to the approach applied in the last section, we obtain 
the scale function using its Wiener-Hopf factorization. It has a form expressed as an infinite sum of exponential 
functions which can be bounded efficiently by finite sums. 

4.1. Meromorphic Levy processes. The following is due to Kuznetsov et al. (2010a), Definition 1. 

Definition 4.1 (spectrally negative meromorphic Levy process). A spectrally negative Levy process X is said to 
belong to the M-class if the following conditions hold. 

(1) The Laplace exponent ip(s) (2.1) has a countable set of real negative poles. 

(2) For every q > 0, the Cramer-Lundberg equation (3.3) has a countable set of real negative roots. 

(3) Let {rjk] k > 1} and k > 1}, respectively, be the sets of the absolute values of the poles and the 
negative roots of (3.3) for fixed q > 0. Then it satisfies the following interlacing conditions: 

• • • < -fik < -£k, q <•••<-%< < -m < < o. 

(4) There exists a > | such that r/^ ~ ck a as k — > oo. 

(5) The Wiener-Hopf factor (3.4) is expressed as convergent infinite products 

The M-class complements the class of Levy processes with phase-type jumps described in the previous section 
because it also contains those of infinite activity. As noted by Corollary 3 of Kuznetsov et al. (2010a), the property 
(3) in Definition 4.1 is equivalent to the condition that the Levy measure has the form (1.1). This can be seen as an 
extension to the hyperexponential case as described in Example 3.1. 

We consider the M-class alternatively to the "hyperexponential fitting". When a Levy measure is completely 
monotone, approximation via hyperexponential distributions is in principle possible. However, as in, for example, 
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Asmussen et al. (2007), special care is needed for the infinitesimal jumps, and one needs to approximate separately 
the process with Levy measure v(—e, 0) for small e > 0. Fitting via the M-class is more tractable in the sense that 
this procedure is not necessary. Although property (4) requires one to choose 77's in a certain way, the approxima- 
tion for the Levy process with a completely monotone density is still effective by choosing the value of c and a 
sufficiently small. For more details, see Kuznetsov et al. (2010a). 

The Wiener-Hopf factor (4. 1) is again a rational function as in (3.5) for the phase-type case. Therefore, this can 
be inverted again by partial fraction decomposition, and we have 



(4.2) 
where 



' { ~K Kq G dx} = A k ^ k>q e-^ x dx, x > 



k=l 



—£ Vq ( S ) 



S=-£fc,9 



6v 



n 



£fc,q 
Vi 



k > 1. 



Vk J " 1 - 



Notice by the interlacing condition that Ak >q > for every k > 1. 

4.2. Scale functions for meromorphic Levy processes. We now obtain the scale function for the M-class. We 
omit the proof because it is similar to the phase-type case; see Appendix A. 1. 

Lemma 4.1. For every q > 0, we have 

oo 

Wq(x) - W Cq (0) = C Kq \l - e-Ko+fc.**"] , x > 

where 



fc=i 



(4.3) 



Cq £,k,qAk,q 



H Cq + £fc,<? 

By (2.1 1) and Lemma 4.1, we have, by taking the limit, 



k > 1. 



(4.4) 



k=l 



The scale function can be therefore obtained by Lemma 4. 1 and (4.4). 



Proposition 4.1. For every q > 0, we have 



(4.5) W^(x)=J2Q, q 



i=l 



i=l 
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By straightforward differentiation, we have, for every q > and x > 0, 

oo 

1=1 

oo 
i=l 

Remark 4.1. From f/ie derivatives above, it can be easily verified that W^' (•) ?j indeed convex. In view of 
Theorem 2. 1, the optimal solution to the classical dividend problem (2.3) is the unique point x such that W^" (x) 
vanishes or 

oo 

i=i 

As an extension to Lemma 3.1, we have the following. 

Lemma 4.2. (1) The following two statements are equivalent: 

(a) a = and u(—co, 0) = oo, 

(b) Efell ^fc.gCfc.g = OO. 

(2) Suppose a > or z/(— oo, 0) < oo. Then, for every q > 0, we have 

(4-6) f = '(f>*^) 

where 



(4.7) := -C q W {q) (0) + (0+) 



771 



A* A 4 

4.3. Approximation of the scale functions via finite sum. The scale function obtained in Proposition 4.1 is an 
infinite sum of exponential functions and in reality its exact value cannot be computed. Here, we obtain bounds 
for W^'(-) and #(«)(•) in terms of finite sums. 

For every m > 1, let 

4™ ) ==l{fc<m} (1-^2) II ST and ^ ): =-?y^' fc > L 



l<i<m,i=£k £ t 



By the interlacing condition, .Aj^ and ^ are all positive and, for every k > 1, 

A^tA fc>9 and t as m -> oo. 
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Now we define candidates for upper and lower bounds of respectively by 



W 



(m) 



i=l 



for every m > 1 and x > 0, where 

m 

i=l 

which vanishes in the limit as m — > oo by (4.4). As candidates for upper and lower bounds of W^ q \ we also define 
W iq,m \x) := e^ x W^\x) and W^ m \x) := e^ x W^\x), x>0. 
The following proposition shows that the scale functions are bounded and approximated by these functions. 



il x _|_ g — (Cq+Cm+l,ij)^ 



Proposition 4.2. For every m > 1 and x > 0, we have 



W_ { l n \x)<W Cq {x)<wT\x) and W^ m) {x) < W {q \x) < W (q ' m) (x) 



(4.8) 

Furthermore, we have 



F (m) 



W £\x) ^ W Cq (x) and W (q ' m \x)^W^(x), 
W^\x)^W Cq {x) and W^ m \x) ^ W^\x), 



uniformly on x £ [0, oo). 



By straightforward calculation, we can bound Zte) in (2.8). Let, for every m > 1, 



Z iq ' m) (x 



:=l + q W {m ' g \y)dy and Z {q ' m \x) := 1 + q / ]£ {m ' 9) (y)dy, x > 0. 
Jo Jo 



Then by Proposition 4.2, we have Z^ m \x) < (x) < Z (q,m > (x) and 

< Z {q ' m \x) - Z (q ' m \x) =q f X (W {q ' m) (y) - W^ m \y)) dy 



q5„ 



o 



x + 



1 



(1-e 



We therefore have the following. 

Corollary 4.1 (Bounds on Z®). We have Z {q,m \x) -»■ Z^(a?) and Z^ m \x) -> zW(i) mm-foo pointwise 
for every x > 0. 
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We now obtain bounds for the derivative. Define, for every x > 0, 

m 

w ( m )( x ) ■— (nh'tf W~ X t cCq® _L H^h 



1=1 



w (m \x) :=w {m) {x) + 

Here notice that 



max (£ k>q e + max (& g e ik < qX ) 

Y<k<m k>m+l 



max (&,e 

k>m+l 



a: ' sm+1,? — x ; 



Proposition 4.3. For every m > 1, we have 

„M (x)<wW(x)< (x), x > 0. 

Furthermore, we have w^ m \x) -4 (x) andw <ym \x) -4 (x) uniformly on x > xq for any xq > 0. 

A stronger result holds when <r > or oo, 0) < oo. Recall in this case that 9 < oo by Lemma 4.2 (2) and 
hence we can define 

e m :=^- k ^^A ( ?>0, m>l, 
which vanishes in the limit as m — )• oo by Lemma 4.2 (2). 
Corollary 4.2. Mien er > or u(—cc, 0) < oo, we have 

u/ m) (x) < W^ Y (x) < w (rn) (x) A w {m) (x) , x > 

w/zere 

:= «/ m )(x) + max (&, g e- & .» a, )<S rn + e- {m + 1 -« a, e ro . 

l<fc<m ,H 

The bounds obtained above on the derivative of the scale function can be used to obtain bounds on the value 
functions and optimal barriers. For simplicity, let w and w be the lower and upper bounds for W^' obtained above 
and let 

u>* = minu>(x) and w* = minu>(x). 

x>0 x>0 

Clearly, w and w are convex because w^ m \ and are for every fixed m > 1. 
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Remark 4.2 (classical dividend problem). Because w and w are convex, we have w* < W^' (a*) < w* and we 
have 

W {q ' m) (x)-5 m (l + e-^) ^ , . ^ W {q ' m \x) 

< u a *{x) < . 

w* w* 

Furthermore, the optimal barrier must lie in the following interval: 

a* G {x > : w(x) < w*} . 

Remark 4.3 (bail-out problem). For the bail-out problem, recall that the optimal barrier is the minimum value of 
a such that G(a) becomes negative (see (2.16)). We can obtain its bounds simply by 

G(a) := [<pZ( q ' m \a) - l]w(a) - <pqW {q ' m) (a) 2 and G(a) := [<pZ {q ' m) (a) - l]w(a) - tpqW< 9 ' m) (a) 2 . 

We can therefore obtain the bounds on the optimal barrier level d*; we have d* < d* < d where 

d* := inf {a > : G(a) < 0} and d* := inf {a > : G(a) < 0} . 

5. Numerical Examples 

We conclude this paper by illustrating numerically the effectiveness of the phase-type fitting approximation for 
a general spectrally negative Levy process. First, we use the classical hyperexponential fitting algorithm for a com- 
pletely monotone density function by Feldmann and Whitt (1998) and approximate, as an example, the scale func- 
tion and solutions for the case with a (Brownian motion plus) compound Poisson with Weibull-distributed jumps. 
Second, we consider, as an example of the meromorphic Levy process, the /3-family introduced by Kuznetsov 
(2009a) and extend the results to the spectrally negative version of the CGMY process. 

5.1. Brownian motion plus compound Poisson process with Weibull-distributed jumps. As noted earlier, 
any spectrally negative Levy process with a completely monotone Levy measure can be approximated arbitrarily 
closely by fitting hyperexponential distributions. Here, we use the fitted data computed by Feldmann and Whitt 
(1998) to approximate the scale function when it is a Brownian motion plus a compound Poisson process with i.i.d. 
Weibull-distributed jumps. Recall that the Weibull distribution with parameters c and a (Weibull(c,a)) is give by 

F(t) = 1 - e~(*/ a ) c , t>0. 

If c < 1, it has long-tails, or e St (l — F(t)) — > oo as t — > oo for any 5 > 0, and has a completely monotone density. 

Feldmann and Whitt (1998) constructed a recursive algorithm to approximate completely monotone densities 
in terms of hyperexponential densities. We use their results and compute the scale functions of spectrally negative 
Levy processes with Weibull-distributed jumps. 
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% oa rji 

1 0.029931 676.178 

2 0.093283 38.7090 

3 0.332195 4.27400 



4 0.476233 0.76100 

5 0.068340 0.24800 

6 0.000018 0.09700 



TABLE 1 . Parameters of the hyperexponential distribution fitted to Weibull(0.6,0.665) (taken 
from Tables 3 of Feldmann and Whitt (1998)). 



Table 1 shows the parameters of the hyperexponential distribution fitted to Weibull(0.6,0.665) when m = 6. As 
can be seen in Figure 4 of Feldmann and Whitt (1998), this fitting is very accurate. We consider the Levy process 
in the form (3.1) where Z is hyperexponential specified in Table 1 as an approximation to Weibull(0.6,0.665). We 
use various values of a with the common values of A = 1, /z = 0.1 and q = 0.03. The roots £., g 's and ( q are 
calculated via the bisection method with error bound 1.0E — 10. 

Figure 1 shows the scale function and its derivative W^' . The optimal barrier levels that minimize W^' 
are given by a* = 0.05, 0.481 and 0.643 for the cases a = 0, 0.2 and 0.4, respectively. The results are consistent 
with Lemma 2.2; with the existence of a diffusion component, the scale function is forced to converge to as a; 
goes to 0. Using these barrier levels, optimal value functions for the classical dividend problem can be computed 
by Theorem 2.1. Figure 2 shows the value functions v n * . Notice that they are monotonically decreasing in a. 
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Figure 1 . Scale functions and their derivatives for the case with Weibull-distributed jumps. 
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o =0 
a =0.2 
o =0.4 



Figure 2. Value function v Wat of the classical dividend problem. 

We now consider the extensions described in the end of Section 2. Here we use the same parameters as in the 
results above. Figure 3 shows the results on the bail-out problem (Avram et al., 2007) when ip = 1.3. It plots G 
in (2.16) as well as the value function v^* in (2.17). Here the optimal barrier level is obtained by computing the 
unique level d* that satisfies G(d*) = 0. Figure 4 shows the results on the extension with terminal values at ruin 
(Loeffen and Renaud, 2010) with the plots of F in (2. 1 8) and the value function %* in (2. 19). We consider the case 
with constant terminal value (i) S = —1 and K = and (ii) 5 = 1 and K = 0. The maximizer of F becomes the 
barrier level b* . Figure 5 shows the value functions v c * on the extension with transaction costs (Loeffen, 2009b) 
when (i) 5 = 0.5 and (ii) 5 = 0.1. In order to obtain the optimal impulse contol (cf , c^), we use the technique 
discussed in Section 4 of Loeffen (2009b). Unlike the other results, the value functions are no longer monotone 
in a unless 5 is sufficiently small. However, as 5 decreases to zero, c\ converges to c* 2 and the value function 
converges to that of the classical model as shown in Figure 2. 

5.2. Numerical results on the /3-class and CGMY process. We now consider, as an example of meromorphic 
Levy processes, the /3-class introduced by Kuznetsov (2009a). The following definition is due to Kuznetsov 
(2009a), Definition 4. 

Definition 5.1. A spectrally negative Levy process is said to be in the /3-class if its Levy measure is in the form 

e a/3x 

(5.1) u(dx) = c _ ^ x l {x<a} dx, xeR, 

for some a > 0, /3 > 0, c>0 and A G (0, 3). It is equivalent to saying that its Laplace exponent is 

j>(z) =p,z+ l 2 a 2 z 2 + 5 l B (a + |, 1 - A) - B(a, 1 - A) 
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n=0 
a =0.2 
a =0.4 



1 -450 
S 



a =0 
a =0.2 
a =0.4 



G(x) 



value function vj* 



Figure 3. Bail-out problem: G in (2.16) and the value function v^* when ip = 1.3. 
Optimal barriers are 0.38, 0.775 and 1.495 for o = 0, 0.2 and 0.4, respectively. 

where B is the beta function B(x, y) := T(x)T(y)/T(x + y). 

The special case a = and (3 = 1 reduces to the class of Lamperti-stable processes, which are obtained by 
the Lamperti transformation (Lamperti (1972)) from the stable processes conditioned to stay positive; see Bertoin 
and Yor (2001) and Caballero et al. (2008) and references therein. For the scale function of a related process, see 
Kyprianou and Rivero (2008). 

It can be also seen that this is a "discrete-version" of the (spectrally negative) CGMY process, whose Levy 
measure is given by 



(5.2) 



v(dx) 



e ax 



dx, x € 



Indeed, if we set c = c/3 A and a = a/3 1 in (5.1), we have 



(5-3) c ( 1 _ e ^)A 1 ^<o> — ► c ]^jX 1 {*<0}' x£R - 

See Asmussen et al. (2007) for approximation of (double-sided) CGMY processes using hyperexponential distri- 
butions. 

We shall use the results in Section 4 to obtain the bounds on the scale functions and the solutions to the classical 
optimal dividend problem. Figure 6 shows the approximation results when q = 0.03, a = 0.2, ft = 0.1, A = 1.5, 
a = 3, P = 1 and c = 0.1 in (5.1). We plot, for m = 15 and m = 150, the upper and lower bounds on the scale 
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o =0 
- o=0.2 
- o =0.4 



a =0 
a =0.2 
a =0.4 



F(x) when (i) S = -1,K = 



a=0 

<j=0.2 

a=0.4 



value function Vb* when (i) 5 



-1,K = 



s 2 
5 



o =0 

o=0.2 

o =0.4 



F(x) when (ii) S = 1,K = 



value function v^* when (ii) S = 1, K = 



Figure 4. Extension with terminal values at ruin: F in (2.18) and the value function iv 
when (i) 3 = -1,K = and (ii) 5 = 1, K = 0. The optimal barriers are (i) 0.0628, 
0.0793 and 0.1384 and (ii) 0.0317, 0.0383 and 0.0955, respectively, for a = 0,0.2 and 
0.4. 



function, its derivative and the function u a * defined in (2.14). As shown in the previous section, the difference 
between the upper and lower bounds indeed converges to zero. 
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Figure 5. Extension with transaction costs: value function when (i) (3 = 0.5 and (ii) 
(3 = 0.1. The optimal impulse controls (cj,c£) are (i) (0,1.173), (0.069,1.527) and 
(0,1.885) and (ii) (0,0.05), (0.222,0.481) and (0.197,0.643) for a = 0,0.2 and 0.4, 
respectively. 

We now take /3 in (5.3) to zero and see how the approximation for the CGMY process works. Here we set a = 3 
and c = 0.1 and use the same values as the above for the other parameters. Figure 7 shows the upper and lower 
bounds of scale function and its derivative for various values of /3. Figure 8 shows the mean value of the upper 
and lower bounds on the function u. Here we can indeed observe the convergence as (3 — > 0. This implies that it 
effectively approximates the scale function and the solution for the CGMY case. 

Appendix A. Proofs 
A.l. Proof of Lemmas 3.1 and 3.2. By (3.6), it is easy to verify that 



8=1 k=l 



x-a 



e ^> qy dy, < a < x, 



and hence, because W^ q ' S C^O, oo), 



i=l k=l 



X, 



(A2) 



d 

Ox 



E x [e- qT0 l 



{ro<oo} I 



x=0+ i=1 
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upper bound 
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upper bound 



Bounds on W^ q '{x) when m = 15 Bounds on W i - q \x) when m = 150 
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Bounds on W^' (x) when m = 15 Bounds on (x) when m = 150 



lower bound 
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lower bound 
upper bound 



Bounds on u(x) when m = 15 Bounds on u(x) when m = 150 

Figure 6. Approximation of the scale function and its derivative for the /3-class. 
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Figure 7. Convergence of scale functions to the CGMY model. 



In fact, different representations of (A. 1) and (A.2) can be pursued. By Theorem 8.1 of Kyprianou (2006) and 
(2-9), 



E* [e-^l {Ta<oo} ] = Z®(x -a)- f W<*\ 

<-><3 



X 



1 + 



W^(y)dy 



W Cq (x 
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Figure 8. Convergence of u to the CGMY model. Here we use M = 150. 



for every < a < x. Its derivative with respect to a becomes 

8 

(A.3) 



.^E* [e-^l {To<oo} ] = -qWb\x - a) + qe^ x ^W Cq (x - a) + j-e^^W^x - a) 



M x - a )w' Cq {x-a). 



In particular, when a = 0, the derivative with respect to x and its limit as x — > are 



(A.4) 

d_ 
dx 



E x [e- qT °l 



{r <oo} 



Q_ 



-C q W iq) (x) + W {qY (x) 



xio+ q 



-C q w^(o) + w^'(o+) 



By matching (A.2) and (A.4), Lemma 3.1 is immediate. 

For the proof of Lemma 3.2, by matching (A. 1) and (A.3) and using Lemma 3.1, we have 



(fc ) (6, g y)* 1 (Ca+ ti 



Qq 



Q q 



i=l fc=l 
n m. 



(fc-1)! 



> 0. 
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Integrating the above and changing variables, we have 



a n m,i ft \ k 

W Cq (x) - W Cq (0) = ^E ^ ' ^ 



(CqHi,q)x 



z k - l e~ z dz, x > 0. 



Lemma 3.2 is now immediate because the integral part is a lower incomplete gamma function. 
A.2. Proof of Proposition 4.2. Notice, for every m > 1, that 

m m 



(A.5) 



and hence by (A.5) 



k=l k=m+l 
m oo 



fe=l 

oo 



k=l 



fc=m+l fc=l fc=l 



fc=l fe=m+l 

Therefore we have the bounds for in (4.8). The bounds for jyw) are immediate by multiplying e 1 ^. Finally, 
the convergence results hold because 5 m — > and e~^ x + e~^" + ^ m+1 < q ^ x and 1 + e _ ^ m + 1 '« :c are bounded uniformly 
in x > 0. 

A.3. Proof of Proposition 4.3. The lower bound is immediate by the fact that < c!-™ < C^ q . For every x > 0, 
we have by (A.5) 



i=l 



w^Hx) 



8=1 



i=m+l 



<w^\x) + max (e M e-^^)^(C iig -C7<^)+ max (^, 5 e^^) £ Q, 9 



Kfc<m 



k>m+l 



i=m+l 



<w( m Hx) + 



max (£ fej? e + max (4 )(? e 

l<k<m k>m+l 



5 m = w^(x) 



which shows the first claim. For the second claim, notice that for the given xq > 0, 



< w {m \x) -w^ m \x) < 



max (£ k , q e~ Ck < qX0 ) + max {i Kq e~^ xa ) 

l<k<m k>m+l 



uniformly on [xq, oo). Because 5 m vanishes as m — > and sup A>0 Ae Xx < oo, the convergence is immediate 



30 M. EGAMI AND K. YAMAZAKI 

A.4. Proof of Corollary 4.2. We have 

m oo 

i=l i=m+l 

m j. oo 

l<k<m 1 — * ' H 1 — * 

i=l ' i=m+l 

< w^ m \x) + max {Ck, q e~^ x )S m + e m e~^ + ^ x 

l<k<m 

where the last inequality holds because 

j. oo j. m j. m 

f E ^A g = d-^Y^^ -jZ^ A ^=^ m>\. 

^ k=m+l ^ fc=l ^ fc=l 

This together with Proposition 4.3 shows the claim. 
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